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ABSTRACT 

The field theoretical renormalization group equations have many common 
features with the equations of dynamical systems. In particular, the manner 
how Callan-Symanzik equation ensures the independence of a theory from its 
subtraction point is reminiscent of self-similarity in autonomous flows towards 
attractors. Motivated by such analogies we propose that besides isolated fixed 
points, the couplings in a renormalizable field theory may also flow towards 
more general, even fractal attractors. This could lead to Big Mess scenarios in 
applications to multiphase systems, from spin-glasses and neural networks to 
fundamental string (M?) theory. We consider various general aspects of such 
chaotic flows. We argue that they pose no obvious contradictions with the known 
properties of effective actions, the existence of dissipative Lyapunov functions, 
and even the strong version of the c-theorem. We also explain the difficulties 
encountered when constructing effective actions with chaotic renormalization 
group flows and observe that they have many common virtues with realistic 
field theory effective actions. We conclude that if chaotic renormalization group 
flows are to be excluded, conceptually novel no-go theorems must be developed. 
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1. Introduction 



The concept of a renormalization group (RG) flow is a basic notion both in 
quantum field theory and string theory [1]-[15]. In the Wilsonian approach, the 
RG equation describes a flow in the space of operators towards a subspacc of 
relevant and marginal operators. This subspace can be viewed as a functional 
attractor for the flow. The behaviour within the functional attractor is then 
determined by the /^-functions for the relevant couplings. Conventionally one 
expects that these couplings flow towards attractors which consist of a finite 
number of isolated fixed points. This ensures that the limit of the RG flow 
yields a definite quantum field theory with destined values for its couplings. 

However, already in [2] Wilson noted that the RG flow of the couplings could 
approach attractors which are more elaborate than plain isolated fixed points. 
In particular, he suggested that the end of a RG coupling flow could be a limit 
cycle. Some evidence in support of this conjecture has been recently provided 
in [16]-[18], where a coupling flow with a (periodic) blow-up is interpreted in 
terms of a limit cycle. 

Here we shall inquire whether a RG flow could indeed tend towards a non- 
trivial attractor, with even a fractal structure corresponding to a chaotic flow 
in the space of couplings. Such chaotic flows could lead to interesting Big Mess 
scenarios in various applications of quantum field and string theories. Indeed, 
we suggest that the commonly accepted dogma that RG flows can only approach 
plain, discreet fixed points is in an apparent contradiction with the existence 
of multiphase systems described by spin glasses and neural networks, which are 
expected to emerge as the IR endpoints of RG flows from simple microscopic 
Hamiltonians. Furthermore, it should be only natural to speculate that the IR 
finality of the full (but yet to be discovered) string theory possesses complex 
multiphase structures [19] with a variety of quantum field theories, strings and 
brane models, emerging at the ends of some as yet unidentified chaotic RG evolu- 
tions. These flows are expected originate from a relatively simple (fundamental 
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microscopic) system, for example from the celebrated M-theory [20, 21]. We 
propose that such chaotic RG scenarios are realizable, very much like condensed 
matter physics with its highly complex long-distance structures emerges from a 
simple microscopic Hamiltonian. Indeed, if complex systems are to appear as 
asymptotic IR limits of some RG flows it should be clear that these flows can 
not be towards a finite set of plain isolated fixed points but rather towards more 
general, even strange attractors with self-similar structures over several orders 
of magnitude. 

We emphasize that we are not considering the possibility of chaos in an un- 
derlying field theory. For example, in the classical Yang-Mills theory chaotic 
behaviour has already been well established [22]. Consequently such chaotic 
behaviour will not be considered here. Obviously, a chaotic RG flow also neces- 
sitates the consideration of field (string) theories with at least three couplings. 

In the present article we shall be interested in the possibility of chaotic RG 
flows in the IR limits of quantum field and string theories. While we are not in 
a position to present concrete examples of theories which exhibit a chaotic RG 
flow, we do have a number of plausibility arguments which support their exis- 
tence. Furthermore, we can understand and explain the difficulties encountered 
in their constructions. At least, our work should motivate the derivation of no- 
go theorems. But at the moment we do not see any immediate contradictions 
between the existence of chaotic RG flows and known properties of field and 
string theories. Indeed, we believe that the potential existence of a chaotic RG 
flow is an important question, and either examples should be searched for or 
then conceptually novel no-go theorems forbidding chaotic RG behaviour should 
be established. 

Exact RG group equations are operative everywhere in the space M. of cou- 
plings {t 1 }. Occasionally, and in particular in the context of string theories, 
these couplings are identified with the moduli space of the theory under con- 
sideration. The Wilsonian approach to the exact RG is based on an effective 
action, 

FA(t\<P) ~ log Z A (t\<p), 
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Z A (t\tp) = j^Dtfiexp ^t'OiWj (with (f>\ aA = <p) (1) 

The integration extends over a functional space A of fields <j>. The effective 
action T depends both on the background fields tp and on a functional form 
which is parameterized by the couplings t l . Consequently it is a section of 
a line bundle over background fields tp and the manifold M. of couplings i*. 
The Wilsonian RG flow (in the sequel we consider mainly flows from the UV 
to the IR) describes the change of T under the change of A, when some of 
the background fields are integrated out. In the simplest case these background 
fields are Fourier components with momenta exceeding some normalization scale 
fi 2 . The boundary dA can have a generic shape. It is parametrized by Whitham 
times T a , and the Wilsonian RG equations can be understood as evolution 
equations in all possible Whitham directions [12, 13]. This leads to a relation 
between the RG flow and the concept of self-similarity (or functional similarity in 
the terminology of [9]) between effective actions evaluated at different Whitham 
times. The present version of renormalizability then makes exact RG flows a 
part of the general theory of dynamical systems [23] . 

The specific feature of a RG flow when viewed as a dynamical system is that 
it involves the effective action. As a section in a line bundle the effective action 
is in general a multivalued function of the background fields and couplings. 
Consequently there is an element of local integrability, that may be eventually 
lost due to global obstructions. Furthermore, since the effective action is defined 
by a (functional) integral it possesses extensive symmetries relating to changes 
of the integration variables (quantum fields) cj>. These symmetries are usually 
expressed in terms of Ward identities (or Picard- Fucks equations) [24]- [26]. For 
an exact RG, the number of independent couplings t must also be large - in 
fact infinite - to ensure that the operators involved indeed form a complete basis 
of functionals. Furthermore, when we ignore the background fields tp, the exact 
RG equations [3] for the effective action acquire a Callan-Symanzik form, which 
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is linear in the derivatives w.r.t. couplings, 



Z(t) = 



dZ(t) 



oz{t) 



dT{t) 



dt l 



^(t) = i i = -ii—, (s = -log M ) 



(2) 



This functional form emerges for a strongly complete basis of operators Oi{<j>) 
(for definitions see [15]) which includes all linearly independent generators. (For 
a weakly complete basis which includes all algebraically independent operators, 
the RG equations contain higher derivatives of Z{t) [15].) 

In the next section we shall first consider certain general aspects of the cou- 
pling flow. We then review some of the assumptions that underlie RG flows, 
i particular the c-theorem(s). In section 3 we first relate a few field theory (3- 
functions to the c-theorems. We then argue by considering the Lorenz model, 
that the properties of these coupling flows are not inconsistent with chaotic RG 
flows. In section 4 we consider limit cycles from the point of view of RG flows, 
and inspect vorticity as a RG scheme independent tool for describing multi- 
coupling flows. In section 5 we study model effective actions as toy models for 
reproducing realistic scaling properties of field theories. In section 6 we explain 
how to construct model effective actions from the /3-function flows. In particular, 
we explain how the construction fails in case of chaotic flows and suggests this 
parallels the problems encountered in constructing actual field theory effective 
actions. This also explains why it is very hard to construct actual field theory 
models with chaotic RG flow. In section 7 we consider possible extensions and 
scenarios for realizing chaotic RG flows, including spectral flow in general and 
in particular in stringy context. We conclude with some suggestions on Early 
Universe models. 
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2. The c-theorems 



The idea of a chaotic attractor is actually not too distant from a quantum 
field theoretical RG equation. This can be seen already by inspecting the func- 
tional form of familiar one-loop single coupling /3-functions. For example, in 
d = 2 + e dimensions the one-loop (i function of the O(N) nonlinear er-model 
has the functional form 

g = eag(l-g) (3) 

which coincides with the form of the Verhulst equation of population growth 
[27] . In its discretized version, this clearly relates to the logistic equation 

x n+1 = cx n (l - x n ) (4) 

which is the classic example of an iterative equation with chaos. 

But for chaos in a continuous RG flow, we need at least three independent 
couplings: In a renormalizable field theory with several couplings t % , the RG 
flow is described by the following autonomous linear system (recall that when 
considering flows towards the IR limit we have s = — In fi which means that the 
ensuing (IR) /3-functions are positive in asymptotically free models) 

-»% = % = ^ < 5 » 

Conventionally, one assumes that in a field theory this IR flow is asymptotically 
approaching an isolated fixed point which is hyperbolic. For classification pur- 
poses we may then evoke the Hartman-Grobman theorem [28] which allows us 
to consider /3-functions which are linear in the couplings, 

(3 l (t) w B'jV + 0(t 2 ) (6) 

(B l j constant). We note that this may not be attainable by conventional changes 
in the renormalization scheme, which are analytic diffeomorphisms on the space 
of couplings of the form 

f -» t\t) = A l fi + 0(t 2 ) (7) 
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where A 1 j is a constant nonsingular matrix with positive eigenvalues. Obviously, 
the first term in (6) describes only integrable, linear flow along the eigendirec- 
tions of the matrix B l j. But it is widely respected that even in the presence of 
apparently very simple non-linear corrections, such as in the case of the Rossler 
system [29] 1 

x = —y — z 
V = x + ay 

z = bx + cz + xz (8) 

the equations (5) can possess a multitude of asymptotic behaviours as the flow- 
time s goes to infinity. Besides isolated fixed points and limit cycles, the tra- 
jectories can also approach chaotic strange attractors with complex geometries 
and fractal dimensions. 

Notice that the equation (5) is renormalization scheme covariant, i.e. form 
invariant under the diffcomorphisms (7). Indeed, under a general coordinate 
transformation t l — > t l (t) the (3 l transform as components of a vector field 

Infinitesimally, for t % = t % + e l (t) we then have 

which is the Lie derivative of (3 along e. Consequently we can interpret the flow 
(5) geometrically in a renormalization scheme independent manner, with /3 l a 
vector field in the tangent bundle of M.. 

In formal quantum field theory investigations one traditionally assumes that 
the flow (5) can only tend towards isolated fixed points. For this, the vector 
field (5 % is subjected to a variety of conditions. A pivotal requirement is that the 
renormalization group flow must be irreversible. According to Zamolodchikov's 

1 Notice that this representation of the Rossler equations differs from the conventional one 
[29] where the last equation reads z = b + z(x + c), by a linear transformation of variables: 
x — * x + ab, y — > y — b, z 2 + 6 so that b = —be and c = c — ab. 
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(weak) c-theorem [30, 31] this irreversibility is ensured by the existence of a 
Liapunov function c(t) which is monotonically decreasing along the RG flow 
towards the 1R 

Occasionally one also assumes that the Liapunov function is positive semidefinite 
c(t) > 0. ft may then be related to the number of degrees of freedom in the 
theory. 

A strong version of the c-theorem states that the vector field (3 l is a gradient, 

/?(*) = -G«(t)^ (10) 

with a symmetric metric Gij(t) — Gji(t). Furthermore, if Gij(t) is positive- 
definite as is usually assumed in the strong c-theorem, the generating function 
<j(t) is a Lyapunov function: 

Furthermore, in the case of two-dimensional field theories it can be argued 
that a(t) is a (peferct) Morse function [32] with its (isolated) critical points 
corresponding to conformal field theories. In that case (11) becomes a gradient 
flow between the critical points of c, i.e. between different conformal field 
theories. 

We note that these statements on the RG flow are renormalization scheme 
independent, and intrinsically geometric. 
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3. Some Examples 



We shall now argue, with examples, that many of the widely accepted prop- 
erties of the RG flow do not exclude strange attractors from appearing in the 
IR limit of the flow. We start by considering examples of field theoretical RG 
coupling flow equations from the perspective of the c-theorems. 

An important example of a RG coupling flow with three independent cou- 
plings (the minimal number required for a chaotic flow) is the U(l) x SU(2) x 
SU (3) standard model in four dimensions. At the two-loop level, and ignor- 
ing the contribution from the Higgs sector, the three gauge couplings t l flow 
according to [33] 

i l = gVdjF (12) 
where the metric g % i has the form 

g V = -SV(b l +J2Vk(t k ) 2 ) (13) 



with 
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where n is the number of generations and F is a (degenerate) Morse function 
for the critical point at t l = 0, 

F(t) = I[(tY-(t 2 ) 4 -(i 3 ) 4 ] (15) 

Here the indefinitcness of F(t) (which is usually excluded by c-theorems due to 
its naive contradiction with unitarity) reflects UV asymptotic freedom of the 
non-abclian components in the model. 
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The first (b l ) term in (13) is the one-loop contribution, the second (b l k) term 
is the two-loop correction. Note that the entire two-loop contribution can be 
viewed as a correction to the one-loop metric. We also note that depending on 
n, the metric can have different signatures. Furthermore, depending on t % the 
signature of the two-loop metric can differ from the signature of the one-loop 
metric. While these observations on the qualitative attributes of the metric as 
such can hardly have much relevance to the Physics of Standard Model they 
are still instructive in revealing the variety of properties that coupling flows in 
quantum field theories share. 

As another example, where already at the one-loop level the signature of 
the metric depends on the relative (small) values of the couplings, we con- 
sider the standard Yukawa coupling between a pseudoscalar meson and a Dirac 
fermion. There are now two couplings, the Yukawa coupling g and the quartic 
self-coupling A of the pseudoscalar. At the one-loop level the flow equations 
are (in d = 4 + e dimensions, with {2-K) d N d the area of the unit sphere in d 
dimensions) [34] 

e 5 o 
9 = 2 9+ 2 d9 + 

A = eA + ^N d (V + ^A.9 2 - 16.g 4 ) + ... (16) 

The origin g = A = is a critical point with its stability depending on the 
sign of e = d — 4. If we introduce the non-degenerate Morse function 

F= 1 -(g 2 + X 2 ) (17) 

these equations can be written in the gradient-flow form 

g = G""d g F = ^(e + 5N d g 2 )d g F 

A = G xx d x F = (e + AN d g 2 - 24N d ^ + '^N d \)d x F (18) 

Obviously the metric can be either positive definite, negative definite or indefi- 
nite depending on the parameters, and the relative strength of the couplings. To 
some extent this can be compensated, by adjusting the relative signs of the two 
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terms in F. But since the metric can also change its signature at non- vanishing, 
even small values of g and A it can not be made positive everywhere. This 
reflects the fixed point structure of the theory on the (g, A) plane. 

As a third example we consider the model which has been studied in [17] 
as a candidate for limit cycle behaviour in the coupling flow. This is the su(2) 
level k = 1 Wess-Zumino model, at the one-loop level its couplings g and h flow 
according to 

9 = -h 2 = -hd g (hg) 

h = -gh = -hd h {hg) (19) 

Consequently in these coordinates the RG equations have the gradient flow 
form but the c-function is not a non-degenerate Morse function. Furthermore, 
unless h > the metric is not positive definite. In the present model the exact 
/3-functions are also known. The ensuing RG flow equations are [35], [17] 

9 = -i^—d g {hg) 
L 9 

h = -J^d h (hg) (20) 
2-5 

In these (isothermal) coordinates we then have a metric tensor which is singular, 
and a c-function which is not a nondegencratc Morse function. Notice in par- 
ticular that the c-function does not receive corrections beyond the one-loop, all 
higher order corrections lead only to a modification of the metric. We note that 
the ensuing coupling flow on the (g, h) plane has a blow-up at finite value of the 
flow parameter s. It has been suggested [16], [17] that this could be interpreted 
as a flow towards a limit cycle (see below). 

From these examples it is clear that the /3-functions that appear in quantum 
field theories are not always of the form suggested by the various versions of 
the c-theorem. In particular, while the equations do admit the gradient flow 
form (10) with (trivially) symmetric metric, the c-functions are not necessarily 
nondegenerate Morse functions nor are the metric tensors necessarily positive 
definite or even regular everywhere in the space of couplings. While in some 
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models such deviations from the c-theorems could be attributed to limitations of 
perturbation theory and removed by higher order corrections, the example (20) 
shows that this is not always the case. Moreover, even though the flows admit 
a gradient representation, they do not in general describe simple, structureless 
laminar flows towards isolated fixed points along gradients of the c- functions: If 
we consider the vector fields which appear on the r.h.s. of the flow equations, 
in all of the three examples, we note that each of the vector fields carries a 
non-trivial vorticity two-form (defined as the exterior differential of the covector 
of the /3-functions) 

Uij = dipj-djpi (21) 

Curiously, for the standard model the vorticity vanishes at the one-loop level 
and appears only at the two- loop level (and beyond). This suggests that higher 
loop corrections have a qualitative effect on the theory. Both in the pseudoscalar 
model and WZW model vorticity is present and regular already at the one-loop 
level. But from (20) we find that the vorticity in the full theory can have a 
singularity at finite coupling 

g(g - 2) - fe 2 
W = = {g -2f 

reflecting the blow-up at finite flow time. We note that since vorticity is a 

closed two-form it can be made constant in a neighborhood of any regular point 

by diffeomorphisms. But if it is non- vanishing in one coordinate system it 

remains non-vaninshing in all coordinate systems. Consequently vorticity is a 

renormalization scheme independent characteristic of the flow. 

Maybe somewhat surprisingly several qualitative aspects of our examples, 

as well as many general assumptions in the c-theorems, can also be realized 

in chaotic systems. As an example, we consider the three dimensional Lorenz 

system [36], 

x = —ax + ay = (3 X 
y = rx — y — xz — (3 y 

z = xy-bz = (3 Z (22) 
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(a, r, b are positive constants). We start by introducing a (Liapunov) c- function, 
an arbitrary positive scmidcfinite function p > in R which we advect along 
the Lorenz flow. This is described by the conservation of the current 

t + v.(^) = 

which reduces to 

Tt = §? +/? ' V/? = -F-ftp = -i* + b + l )p < (23) 

Consequently any positive function on R 3 which is advected by the Lorenz flow 
satisfies the irreversibility requirement of the c-theorem. 

Consistent with the c-theorems, the Lorenz system can also be presented in 
the form of a gradient flow (10) with a symmetric metric. For example, if we 
introduce a non-degenerate Morse function for the critical point at the origin 

F(x,y,z) = x 2 + a(y 2 + z 2 ) (24) 



we can write the Lorenz system as 

x = g^dxF = a V —^d x F 
2x 

y = g 22 d y F = rX ~ y ~ XZ d y F 
y J y 2ay 11 

i = g ^ dz F = X -V-^d z F (25) 
Aaz 

As in the previous three field theory examples, the positive definiteness of the 
metric tensor g a b depends on the relative values of x, y and z. Like the metric 
in (20), this metric is also non-singular except for the critical points of the flow 
(and the x, y, z — lines) while F is both a Morse function for the critical point 
at origin, and a global Liapunov function when r < 1 since 



dF n , 

— = -2a ( 
dt y 



+ -(l-r)(3 + r)y 2 + 2abz 2 ) < {x,y,z^0) 



Consequently, at this level of analysis we do not see much difference in the 
qualitative properties of the (chaotic) Lorenz model and the coupling flows in 
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our three field theory examples. In particular, nothing appears to prevent chaos 
from occuring in RG flows as well. 
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4. Limit Cycles and Vorticity 



The existence of a positive definite metric with a c-function which decreases 
along the orbits can also be satisfied by flows without simple fixed points. For 
example, if we take 



the entire cycle t = ae 10 forms an attractor. In order to induce a motion along 
this cycle we then consider (rj > 0) 

x = —r\y = —(3 X 



This flow is consistent with the strong c-theorem, with a positive-definite metric 



and a (multivalued) Liapunov c-function which is monotonically decreasing 
along the flow, 

c(x,y) — arctan(y/x) 
The ensuing flow has constant vorticity, 

cj = codify = 2r] 

Such flows may actually be realistic in models where the coupling (moduli) space 
has nontrivial topology, with non- vanishing •ki{M.) (cf. [37]). Indeed, since the 
c-function is highly nonlinear the flow could be viewed as a non-perturbative 
one. We are not aware of any apparent reason why this kind of flow should in 
general be excluded by the c-theorem. 

There is a simple generalization of (26) to non-constant vorticity, 



c(t,t) = {tt-a 2 ) 2 



y = r]x 



/3 y 



(26) 



= ^ (x 2 + y 2 ) 




y = x 1 



(27) 
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with k an integer. This flow is of the form of strong c-theorem, with 

c(x,y) = y Y - k -x Y - k 
The corresponding metric tensor is 

a 13 = —^y) 

For k = 2n this metric is positive on the (x, y) plane, and the c-function has the 
form of a (degenerate) Morse function as expected by the strong c-theorem. But 
for k = 2n+ 1 i.e. odd, the metric is in general not positive. It is instructive to 
consider in more detail the reasons for the failure of strong c-theorem when k 
odd. For this we note that the flow possesses a conserved quantity, x k+1 +y k+1 = 
const. This suggests us to change variables 

x = i?cos 1 /(™ +1 ) e 

y = i?sin 1/( ™ +1) 6 (28) 

so that 



R = 0, 
dc{9) 



(29) 



86 

where c(0) is R- independent. Consequently, in these coordinates it appears that 
the flow is consistent with the strong version of the c-theorem with a positive 
definite metric 

dR 2 + R 2 d8 2 

But when we transform back to the cartesian x, y coordinates, 

, . z dXdY dXdY\ 

(30) 

We find that due to the fractional powers of trigonometric functions in (28) the 
ensuing metric is not positive definite, for example when k = 1 we have 

qXX ( R S m0 \ 2 | R 2 (31) 





2\/cos 9 ) \ 2 I 4cos0 



16 



which is not positive definite. 

Notice that besides a gradient flow, the flows (26), (27) can also be inter- 
preted in terms of a symplectic (Hamiltonian) flow of the form 



with a closed symplectic two- form ut = uj^dt 1 A dP and a Hamiltonian H. In 
particular, (26) is the harmonic oscillator. 

We note that the presence of such limit cycles in RG flows is quite impor- 
tant from the point of view of (possible) chaotic flow. Obviously, a limit cycle 
behaviour is much easier to analyse than a chaotic behaviour. But in addition, 
limit cycles can also provide a period doubling (Fcigcnbaum) route to chaos 
(for flows with at least three couplings) [38]. Indeed, if for some value of ex- 
ternal control parameters a system exhibits a stable limit cycle, its stability 
can be lost by period doubling when the control parameter changes. When 
this happens, the attractive cycle becomes repelling and instead there is a new 
attractive limit cycle which exhibits period doubling and links around the pre- 
viously stable cycle. When the control parameter is varied further, there will 
be additional period doublings and eventually these bifurcate into an infinitely 
long limit cycle, with a fractal structure. The Rossler system (8) is the simplest 
three dimensional example which exhibits this Fcigenbaum route to chaos by 
consequtive period-doublings in its limit cycles, while in the Lorenz system (22) 
the transition to chaos is by intermittency [27]. 

Unfortunately, it appears to be quite difficult to find RG limit cycles. While 
periodic dependency on the coupling constants in the /3-functions has been well 
established, for example in the context of topological charge rcnormalization 
[39], this is not sufficient for obtaining a limit cycle. Until now, only one example 
of limit cycle behaviour in RG equations has been constructed [16]-[18], but 
unfortunately this example can not be considered fully satisfactory. Essentially, 
the flow is of the form 




(32) 



g = h 2 +g 2 , 
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h = 0, 

which exhibits a blow-up, rather than smooth cyclic behaviour. We note that for 
small values of g these equations are quite similar to the conventional asymp- 
totically free RG equations (with h = 0), the only difference is that whenever 
the second coupling constant h ^ the flow becomes accelerated. The ensuing 
flow 

g{n) = h cot (h log p) 
is formally periodic in the sense that 

g{e*' h p) = g{ii) 

but a discontinuous jump from g = +00 to g = —00 is necessary in order to 
close the cycle. Unfortunately concrete examples with continuous limit cycle 
behaviour in realistic field theory models have not yet been found. 

The previous examples underline the importance of developing general clas- 
sification schemes for /3-function flows in multidimensional cases. As we have 
proposed, a local approach could be based on the Hartman-Grobman theorem 
[28] on hyperbolic fixed points. This suggests that for classification purposes 
we consider flow equations with /3-functions which are at most bilinear in the 
couplings, the bilinearities representing cither corrections to, or deviations from 
hyperbolic behaviour. This subclass of flows is particularly interesting in three 
dimensions, since bilinear nonlinearities are sufficient for a chaotic flow to occur. 
For this we consider three dimensional flows of the form 

x* = f3 l (i = 1,2,3) (33) 

The three dimensional (co)vector can be presented using a Glebsch decom- 
position 

fa = fidip + dij 
with three functions jj,, p, 7. When \x = the vorticity 

= eijkdjPk = Eijkdj/idkp 
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vanishes, limit cycles are absent and the ensuing advection of the couplings 
is laminar, non-chaotic gradient flow. But whenever u is non-vanishing the 
vorticity is non- vanishing and either constant or linear in the couplings. We now 
argue that both for constant and linear vorticity, the advection can be chaotic. 
(For a flow with vanishing vorticity, chaotic advection is hardly possible.) Recall 
that since vorticity is a closed two-form, we can introduce a diffcomorphism 
which maps it into a constant in a neighborhood of a regular point. But it can 
not be made to vanish in that neighborhood by diffeomorphisms. 

Consider first the Lorenz equations (22). The (3- functions can be Glebsch- 
decomposed according to 

Pi = xy'di ( i± ^- 1 ) + 9i [axy \{ox 2 + y 2 + bz 2 )^ 

and for the vorticity we find 

lo = —zdx A dy + (2x — a)dy A dz — ydz A dx (34) 

Notice that on the surface x 2 + y 2 + z 2 = const this involves a term which rep- 
resents H 2 (S 2 ), the (unique) volume element of S 2 . Consequently the vorticity 
(34) of the Lorenz system is a representative of the monopole bundle in R 3 . 

For the Rossler equation (8) we find the vorticity (here we use the original 
form of these equations, see footnote (1) in connection of equation (8).) 

u) = -2dx Ady+(1 + z)dx A dz = -2dx A dy + ^dx A d(l + zf (35) 

Consequently the vorticity can be made constant by a quadratic diffeomorphism, 
but at the expense of loosing the quadratic nature of the equations. 

A chaotic system with simultaneously quadratic nonlinearities and constant 
vorticity can also be constructed, for example by combining the behaviour of 
the /3-functions in (19) and (26) into the following three dimensional flow 

1 1 2 

Px = -^z + y-l+px+-y 

P y = -x + xy - y - 1 

P x = -^x + az (36) 
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This can be shown to be chaotic for appropriate values of the parameters p, a 
[40] . Vorticity is linear 

uo = —2dx Ady + dz A dx 

In D dimensions we can introduce a generalization of the Glebsch decompo- 
sition 

N 

k=i 

Where 2N = D for D even, and 2N = D — 1 for D odd, and 7 is present only 
for odd D. Presumably a necessary condition for chaotic advection in D > 3 is 
that the ensuing vorticity 

Wij = ^{diUkdjPk - djpkdiPk) 
fe=i 

is non-vanishing. Since a non-trivial vorticity appears to be generic for (3- 
functions in quantum field theories, additional restrictions are then needed to 
exclude chaotic advection of the couplings. 
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5. RG Flows And Model Effective Actions 



= -26, (37) 



It appears that nontrivial vorticity is generic for the coupling flows, to the 
extent that even the one-coupling /^-functions can be related to flows with vor- 
ticity. For this, we consider the familiar one-loop formula 

d_ ( l_ 

ds \g 2 , 

which appears e.g. at one-loop four dimensional Yang-Mills theories, cf. (12), 
(13). We then introduce a Coleman- Weinberg type model effective action (for 
strong fields). This is an ordinary function of a single real variable F, 

T(g,h\F) = ±F 2 + hF 2 log F, (38) 

It turns out that many properties of RG flows can be understood by inspecting 
the scaling properties of such ordinary functions. Indeed, since the /3-functions 
in (5) are independent of spacetime coordinates, we can expect that the essential 
aspects of RG flows are independent of spacetime coordinates and can be studied 
in terns of such ordinary functions (i.e. model effective actions) in lieu of the 
actual quantum effective actions. This is certainly the case in theories where 
we can have an effective potential, such as Higgs models. 

When we introduce the scaling F X 2 F and T — > \ A T in (38), we get 

ds (s 2 ) g 3 

h = (39) 

The solution 

h = const = —2b 
g~ 2 = 2b log a 

of these equations leads to the /3-function (37). Nevertheless, the system (39) 
has a non-trivial vorticity 

w = dhhdg- 2 ± (40) 
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More generally, if we consider model effective actions of the form 

T{h\F) = Y J hkF 2 \og k F (41) 

fe=0 

the scaling F — > \ 2 F leads to the evolution equations 

h k = {k + l)h k+u (42) 

or, for g k = h k /h Q , 

g k = (k + l).gfe+i - gigk- (43) 

If all h k with k > N vanish, (42) is solved by a polynomial Pn(s) of degree 
N for ho, and its fc-th derivative for h k . Thus solutions g k of (43) are rational 
functions of the form g k = d k Pn (s) / Pn (s) ■ The scaling properties of (41) and 
(42) can then be employed to derive realistic multiloop /3-functions of the form 

m = E b ^ 2j+1 

3 

generalizing our one- loop result (38). For this, we substitute for h k in (41) 

1 

7 



ha(g 2 ) - .,■ 



^J-*.— nS^-'. 



9° 9 



h 2 (g 2 ) =h 1 = 2gf3(g)h[(g 2 ) = -4^(j - 1^9^+^, 



(44) 



If & = 0, then 

/ii = -2(6! +b 2 g 2 + ...), 
h k = -2 k b\- x b 2 k\g 2k + ... for k > 1. (45) 

The scaling properties of the model effective actions (38), (41) lead to the 
Callan-Symanzik type equations 

« + E^)« + W = (46) 
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where the model background field ip in (46) is related to F in (41) by F = 
e~ v . The equation (46) can then provide relations between field theoretical 
RG equations and dynamical systems. In particular, we shall now propose that 
issues concerning chaoticity and integrability of the coupling flow equations can 
be directly related to the construction of solutions to the equation (46). 
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6. Model Effective Actions And Chaos in RG Flows 



In quantum field theories, it is known that the effective actions are usually 
highly complicated, multivalued functions of the couplings. Indeed, in general 
the effective actions are not ordinary function(al)s but highly nontrivial sections 
of line bundles which are defined over fibrations of the background fields <p over 
the spaces M. of couplings ij. We shall now propose that these complexities in 
their functional form could be viewed as an indication of chaotic behaviour in 
the ensuing flows of the couplings. This turns out to be an issue that can be 
addressed at the level of model effective actions i.e. ordinary functions which 
model the scaling properties of the actual field theoretical effective actions. 
These model effective actions can often be constructed by explicit integrations 
of the ensuing RG flow equations: 

Consider the model RG equation (46), which we now solve using the method 
of characteristics. For this, we first need a solution to the extended system (5), 

h=(3i(t), 

(p = 1 (47) 

These solutions yield the flows ti(s,cti-i), with dj_i as the integrations con- 
stants. We then invert these relations, to express s and the integration con- 
stants a,_i as functions of the tj and ip. The general solution to the model RG 
equation (46) is now obtained by introducing an arbitrary function /[dj_i] of 
di-i{t\(p) and setting 

F(t\<p) = f[a i - 1 (t\<p)]+k8(t\<p). (48) 

Obviously, complex behaviour such as bifurcations in the original equations (47) 
relate to singularities in these functional inversions. 

Clearly, if the equation (5) describes flow towards a strange attractor i.e. 
the flow is chaotic, the solution (48) will also lead to a function which reflects 
the structure of the attractor. Asymptotically, the model effective action (48) 
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then flows towards a generalized function or rather a (Lebesque) measure with 
support on the strange attactor. Obviously this means that we need to extend 
the concept of model effective actions to include measures with support on 
fractal structures. In particular, when the flow (47) is chaotic the construction 
of its solutions can only be symbolic, which then translates to the impossibility 
to construct the model effective action (48) by quadratures. Presumably, this 
could be developed into a criteria for identifying actual field theory models where 
chaotic RG flows are present. 

For this, we first consider explicit constructions of model effective actions 
for the flows that we have analysed previously. We start by noting that the 
last equation in (47) is clearly oversimplified and easily solved, ip = s + a s , so 
that ip will actually enters (48) only through a single integration constant a s . 
Moreover, it will always appear in the combination a s = ip — s(t), where s(t) is 
obtained, together with the integration constants aj_i, as the functional inverse 
of the flows, derived from (5). For the same reason the last term at the r.h.s. 
of (48) does not depend on tp and can be written as ks(t). 

We first consider the flow of a single coupling, 

t - /?(*) (49) 

The generic solution to (46) is 

F{t\<P) = fW) -<p] = f [7We-"] , (50) 

where f[x] = f[e x ] is arbitrary function of a single variable, and 7 ; (i)/7(i) = 
l//3(i), i.e. 

7 ( *) = ^=exp(/^) (51) 
Similarly, if the i-th (3- function depends on only, (3i(t) — 0i(ti), 

m<P) = fbi(U)e- v ], (52) 
where f[x] is arbitrary function of its variables and 

^^{Pm) (53) 
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We then consider 



ti -A-ij tj j 



(54) 



with i-indepcndent matrix Ay. According to the Hartman-Grobman theorem 
[28], a generic flow is homeomorphic to this near its hyperbolic fixed points. 
But we note that this can also include models with limit cycles, such as (26). 
We diagonalize Ay so that the system is transformed into & = with 



(UAU-% = XAj 

The generic solution to (46) then becomes 

a/A; 



HtW) = f 



f''"* i 

-Tn-, <p- t- log a 



-\np 



(55) 



with arbitrary functions and /[x,]. 

Finally, we consider the most general linear flow 

t{ A^ tj -\- Bi , 

with t-independent A^ and B t . We get 



HAv) = f 



(56) 



(57) 



Note that we have defined the arguments of f[xi] so that they make sense even 
if A^ is degenerate i.e. some of the eigenvalues Ai vanish. 

A particular example of the previous construction is given by the limit cycle 
flow (26), 



x = -rjy, 
y = i]x 



(58) 
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we have 



-asmr/s, y — a cos rjs, and the generic solution is 



F{x,y\<p) = f V x2 + V 2 , V<P + arctan(a;/y) = 
= / [x cos r]<p — y sin r](p, x sin r)<p + y cos r)ip] , 



(59) 



where f[x,y] and /[a;, y] are arbitrary functions of two variables. 

As a further, nonlinear example we consider the model effective action for 
the /3-functions (19), (20). For this we consider (19), in the form 



x = xy, 



y = x 



(60) 



Indeed, this is the system considered in [16, 17], with 

Q 



x{s) = — 

sin Ljs 

The generic solution of (46) is 



& y(s) = Q cot Qs 



Qs'mQif - ycosQf 



(61) 



where f[x,y] and f[x,y] are arbitrary functions of two variables, and 



Q(x,y) = yjx 2 - y 2 

It is clear that since the method of characteristics is based on solving the 
system (5), in the case of chaotic flows the method can not be explicitly im- 
plemented, even in principle. Consequently model effective actions for chaotic 
flows can not be explicitely constructed, which may also be an explanation why 
concrete examples are not known in the literature. Indeed, this is an obvious 
conceptual issue that explains why it is very difficult (may be even impossible?) 
to construct even simplistic model effective actions with chaotic renormaliza- 
tion group flow in the first place. To exemplify the problems encountered, we 
shall now consider the construction of the model effective action for the Lorenz 
equations (22), by employing a perturbative expansion. A natural perturbative 
parameter is r, which corresponds to the ratio of Rayleigh number to its critical 
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value in the underlying hydro-dynamical model. This suggests that we search for 
a perturbative construction of the model effective action by first expanding (22) 
around the (chaotic) large-r limit. (Recall that with the canonical values a = 10 
and b = 8/3 the model becomes chaotic only when r > 24.74 ...) The expansion 
parameter we use is e = 1/y/r, and setting x — > ex, y — > e 2 ay, z — > <r(e 2 z — 1) 
and t — > t/e we get for the Lorenz equations 

x = y — eax 
y = —xz — ey 
= xy-eb(z + a) (62) 

x + exi + e 2 x 2 + ... 

yo + m + £ 2 y2 + ••• 

z + ezi + e 2 z 2 + ... 

xo = yo 
yo = "X z 

zq = x a y 

These equations can be integrated by quadratures, in terms of the Jacobi elliptic 
functions. When we substitute the solution to the higher order equations for 
{x n , y n ,z n ) (n> 1) we find at each order a set of equations which are linear in 
their unknown variables. For example at (9(e) we get the linear equations 

Xi =2/1- crx 
yi = -x Q zi - xiz - y 
zi = x yi + xiy - b(z + a) 

and so forth at higher orders in e. Consequently these higher order equations can 
in principle also be solved by quadratures. Thus the Lorenz equations can be 



We then expand 



x 

y 

z 



To the leading order we get 
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solved by quadratures, at least formally in a perturbation expansion to an arbi- 
trary order in e. The ensuing model effective action can then also be constructed 
by employing the method of characteristics, order by order in a perturbative ex- 
pansion in e. However, it turns out that this perturbative solution of the Lorenz 
equations is at most asymptotic. It does not converge towards a solution of the 
original Lorenz system. Instead, numerical investigations indicate that it leads 
to an diverging asymptotic expansion which describes the chaotic solutions of 
the Lorenz equations only for (relatively) small values of the flowtime s. 

Clearly, it must be a general feature of chaotic systems that model effective 
actions can not be constructed by quadratures. Not even in principle, as this 
would amount to solving the original chaotic equations. For the same reason 
any perturbative approach can only lead to an asymptotic expansion, which ap- 
proximates the exact model effective action at most during a limited period of 
the flow. Indeed, since a chaotic flow approaches an attractor which is a fractal, 
the ensuing model effective actions must also approach generalized functions 
(measures) with support on a fractal set of points. We note that this is in a 
very curious resemblance with the familiar complex behaviour of actual effec- 
tive field theory actions. These effective actions are usually highly complicated 
and multivalued function(al)s (rather section(al)s) for which any perturbative 
expansion yields at most an asymptotic series. 
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7. Further Developments 



It was proposed already in [2] that (at least) cyclic behaviour of asymptotic 
RG flows can not be excluded. As we have argued in the present article even 
chaotic flow does not appear to conflict (most of) the assumptions in c-theorems. 
In fact, cyclicity of the flow seems to be quite natural in many field theory prob- 
lems with spectral flow, with an adiabatic evolution of the energy eigenvalues. 
In terms of spectral flow, a cycle can arise whenever the energy spectrum is 
mapped onto itself by an adiabatic process, but with a rearrangement: The 
level E n becomes shifted into some level E n -k- From the point of view of an 
effective action, this rearrangement of the energy levels is clearly a cycle. 

One example of nontrivial spectral flow can be developed by considering a 
stringy spectrum of the form 

E n =nm 2 (g)+a(g) (63) 

A cycle arises whenever 

m 2 {g') = m 2 (g) and a(g') = a(g) - km 2 (g), k e Z. (64) 

The intercept a(g) is now identified as the dissipative Lyapunov c-function. 
Note that in simple string models the intercept is exactly (proportional to) the 
central charge of the underlying 2d conformal field theory (effective space-time 
dimension). Clearly, this example is very much in the spirit of the original 
introduction of the c-theorem in the context of conformal field theories [30] . 

In our examples, we have inspected the c-theorem in RG flows without 
background fields. To some extent this can be justified, by arguing that the 
background fields can be amalgamated with couplings. However, there are also 
differences. For this, we consider a hypothetical field theory model ( for example 
a Higgs model), with background fields which describe a vacuum state of the 
effective theory. If the parameters in the effective action change e.g. as a conse- 
quence of coupling flow, then the functional form of the effective action will also 
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change. But as the effective action changes, so does the vacuum. Moreover, 
in the case of first-order phase transitions there are discontinuous changes in 
the ensuing order parameters i.e. background fields. Even if only second-order 
phase transitions occur in the course of RG evolution and the location of vac- 
uum is changing continuously, its evolution is still important for determining the 
effective action. The physical partition function which emerges after accounting 
for the change in the vacuum obeys RG equations somewhat different from the 
original theory. For this, we consider a model effective action 

k 

which flows as 

Fif) = B {<P) = ( 65 ) 

k 

Define a vacuum order parameter tpo as an extremum of F{(p): 

F'{<Po) = 0. (66) 
Often we are interested in effective potential 

f„( 9 )E%„^)-jy, (67) 

which describes fluctuations around the vacuum. Then, since 

= -B(<p )/F"(<p ) 



we get 



and 



F {ip) = B (<p), 



B {ip) = Bfa + v) - B(v ) - +<p)~ ^'M] (68) 

Now, even if B(<p) is subject to a gradient flow according to the strong c-theorem, 
there is no apparent reason why this should be the case with B (ip). 
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Once we account for the RG flow of the vacuum, new phenomena become 
possible. For example, consider the potential 

V{<f>) = + ^ 

The vacuum </> = ^Jm 2 / X changes smoothly when to 2 and A are renormalizcd 
according to conventional RG equations. Now suppose the field cf) is further 
coupled to another field \i f° r example by 

X 4 cos(w<£) 

Then the effective coupling 34 in the x 4 vertex is actually equal to 

.94 = cos(w</)o) = cos(o> \J 'to 2 / 'A) 

and, for sufficiently large u>, this (74 can be oscillating along the RG flow with 
ensuing changes for to 2 and A. This model could then indicate how to construct 
realistic field theory models with RG limit cycles. 
An alternative is to consider 

V(4>) = g\ cosu)\(j) + 52 cosw20 (69) 

with u>i/u)2 irrational. The potential has infinitely many irregular local minima, 
which change when gi and gi are smoothly evolving along any RG flow (which 
may be subject to the strong c-theorem). The lowest energy minimum is a 
complex function of these couplings, and the ensuing v.e.v. <po can at least a 
priori exhibit chaotic (irregular) behaviour. Such irregular behaviour of (j) can 
then give rise to similarly irregular behaviour of some of the couplings. 

Similar phenomena can also be realized in the context of finite tempera- 
ture field theories, with various applications including in particular Very Early 
Universe and Cosmology. In the presence of a finite temperature, the role of 
flow-time s is taken by the inverse temperature s = l/T. As temperature de- 
creases, the number of states that contribute to the partition function decreases 
- the flow is contracting, irreversible. Again, we can illustrate the phenomena 



32 



with a model partition function: Wc select a somewhat unconventional set of 
harmonic oscillators, with potentials 

The ensuing spectrum consists of oscillators +uiiki with integer non-negative 
fcj, with the ground state energy absorbed into the intercept on = a, + The 
partition function is 

Z{s) = exp(^(.s)) = ]T ( ? °) 

If aj ^> Wj, this can be approximated by 

1 ^ - e -«( w )» 

*w«) = -E— ^— (71) 

In the presence of only two frequencies uj\ <C 0J2, with ai » a 2 , we then have 
sZ app (s) = + , (72) 

Ul\ 0l>2 

with the first term dominating at s -C 1/ai and the second one at 1/ai «s« 
l/«2- For infinitely many chaotically distributed LUi and ai, the behaviour of 
Z(s) also exhibits irregular (chaotic) features. As an extreme, one can consider 
a somewhat peculiar distribution of oscillators with to n ~ y/n and a B = a log n. 
Then 

sZ app (s) = C Q + as) (73) 

is the Riemann ^-function, which is suspected to have a relation with a chaotic 
dynamical system. 
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8. Conclusions 



In conclusion, we feel that there is a clear need to investigate the relations 
between RG flows from the perspective of chaotic dynamical systems. We have 
analyzed the restrictions that can be imposed on the RG flows by employing 
various aspects of c-theorems and local intcgrability conditions, stemming from 
the fact that RG flows reflect the scaling properties of effective actions. These 
effective actions are in general not functions but sections in a line bundle over 
background fields and couplings. In fact, we have argued that the involved 
structure of field theory effective actions, in particular their multivalucdncss 
in the couplings, can be naturally understood in terms of the impossibility to 
cxplicitely construct model effective actions for chaotic dynamical systems. In 
particular, both admit perturbative expansions which fail to converge except in 
an asymptotic sense. 

We have also proposed that many of the familiar properties imposed on 
RG flows arc not in any kind of apparent contradiction with the existence of 
chaotic behaviour. In fact, several RG flows do reflect features such as nontrivial 
vorticity which are necessary for a chaotic flow to occur. Furthermore, we have 
suggested that a self-similar, chaotic RG flow in the IR limit could actually 
be desirable in many applications, e.g. to spin-glasses and neural networks. 
Perhaps even more importantly to the Early Universe Cosmology, and maybe 
even M-theory with the possibility that various field theory, string and brane 
models emerge at the end of its chaotic RG trajectories. In a sense, we are then 
proposing that chaotic RG orbits are as natural as the emergence of condensed 
matter physics with its highly elaborate and largely self-similar structure which 
emanates from the relatively simple microscopic quantum electrodynamics. 

While we feel that the idea of self-similar chaotic orbits in field theory RG 
flows is a very natural one, we can not exclude their absence. But for this, novel 
no-go theorems are needed. Such theorems should most likely be based on 
conceptually new physical principles. Whatever the necessary restrictions are, 
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they arc bound to shed important light to studies of hidden integrable structures 
of effective actions (with complex relations between different RG flows taken 
into account). They will also lead to a deeper understanding of singularities of 
generalized r-functions, and in a wide sense to the general structure of phase 
transitions with applications ranging from condensed matter to early universe 
and fundamental string theories. 
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